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We report on recent studies of the spin-half Heisenberg and the Hubbard model on 
the sawtooth chain. For both models we construct a class of exact eigenstates which 
are localized due to the frustrating geometry of the lattice for a certain relation of the 
exchange (hopping) integrals. Although these eigenstates differ in details for the two 
models because of the different statistics, they share some characteristic features. The 
localized eigenstates are highly degenerate and become ground states in high magnetic 
fields (Heisenberg model) or at certain electron fillings (Hubbard model), respectively. 
They may dominate the low-temperature thermodynamics and lead to an extra low- 
tcmpcrature maximum in the specific heat. The ground-state degeneracy can be calcu- 
lated exactly by a mapping of the manifold of localized ground states onto a classical 
hard-dimer problem, and explicit expressions for thermodynamic quantities can be de- 
rived which are valid at low temperatures near the saturation field for the Heisenberg 
model or around a certain value of the chemical potential for the Hubbard model, re- 
spectively. 

Keywords: frustration, Heisenberg model, Hubbard model, localized eigenstates 

1. Introduction 

Frustrated lattices play an important role in the search for exotic quantum states 
of condensed matter. The term 'frustration' was introduced in physics in the 1970s 
by Toulouse [1] in the context of spin glasses [2] and describes a situation where 
exchange interactions are in competition with each other. The studies on spin glasses 
have demonstrated that frustration may have an enormous influence on ground-state 
and thermodynamic properties of spin systems [2] . 

In the 1970s Anderson and Fazekas [3] first considered the quantum spin-1/2 
Heisenberg antiferromagnet on the geometrically frustrated triangular lattice and 
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proposed a liquid-like ground state without magnetic long-range order. Although 
later on it was found that the spin- 1/2 Heisenberg antiferromagnet on the triangular 
lattice possesses semi-classical three-sublattice Neel order (see, e.g., Refs. [4,5] for 
recent reviews), Anderson's suggestion was the starting point to search for exotic 
quantum ground states in frustrated spin systems. 

The recent progress in synthesizing frustrated magnetic materials with strong 
quantum fluctuations [6] and the rich behavior of such magnetic systems have stim- 
ulated an enormous interest in frustrated quantum magnets, see, e.g., Refs. [7-11]. 
There are many compounds which correspond to quantum antiferromagnetic 
Heisenberg models with frustrated spin interactions. We mention as examples the 
frustrated spin-1/2 Ji- J 2 chains (Rb 2 Cu2Mo 3 Oi2, LiCuV0 4 , Li 2 ZrCu0 4 ) [12] and 
the kagome lattice (ZnCu3(OH)6Cl2) [13]. There are also compounds which corre- 
spond to electronic (Hubbard, t — J , periodic Anderson) models on geometrically 
frustrated lattices. We mention as examples cobaltates [14], CeRhjB^ [15], as well 
as artificial crystals from quantum dots [16]. 

In this paper we will focus on a special property of the Heisenberg and the 
Hubbard model on a particular geometrically frustrated lattice (the sawtooth chain, 
see Fig. [IJ , namely the existence of localized eigenstates (on a perfect lattice) and 
their relevance for the low-temperature physics of those correlated systems. Note, 
however, that arguments and calculations presented in this paper can in principle be 
applied to wide class of frustrated lattices, see the discussion below and Refs. [17-28]. 

In general, for perfect lattices an elementary excitation as a non-interacting 
quasiparticle is spread over the entire lattice. For example, for a simple hypercu- 
bic lattice a magnon or electron wave function is extended over all lattice sites 
due to a hopping term in the Hamiltonian. However, for some lattice geometries a 
wave function of an elementary excitation in a quantum system may have ampli- 
tudes which are non-zero only in a restricted area owing to destructive quantum 
interference. We call such excitations localized excitations (for example, localized 
magnons [17-19] or localized electron states [25-30]). Due to the local character of 
these excitations exact many-particle eigenstates of the Hamiltonian can be built 
by n independent localized excitations (i.e. they have a sufficiently large separation 
between each other) even in the presence of interactions. The number n of local- 
ized excitations cannot exceed a certain maximal value n max which depends on the 
specific lattice under consideration, where n max is proportional to the system size 
N [18,28]. If the localized excitation is the lowest-energy eigenstate of the Hamilto- 
nian in the one-particle subspace one may expect that a state with n independent 
(isolated) localized excitations is the lowest-energy eigenstate of the Hamiltonian in 
the corresponding n-particle subspace [17,18,28,31] provided there is no attractive 
interaction. The localized eigenstates may become ground states in high magnetic 
fields (Heisenberg model) or at certain electron fillings (Hubbard model), respec- 
tively. Therefore they may substantially contribute to or even completely dominate 
the low-temperature thermodynamic properties of the system. 

In the present paper we discuss the effect of localized elementary excitations on 
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Fig. 1. (Color online) Upper part: the sawtooth chain. Filled circles indicate the lattice sites, lines 
indicate the exchange/hopping bonds. Two trapping cells occupied by localized magnons/electrons 
are indicated by bold lines. Note that for the sawtooth chain one has two kinds of bonds of different 
strength. The lower part of the figure indicates the corresponding hard-dimer model (two hard 
dimers on a linear chain). 



the low-temperature thermodynamics focusing on the quantum Heisenberg antifer- 
romagnet and the Hubbard model on the sawtooth chain. We follow the lines which 
have been developed in a series of papers on localized eigenstates for the Heisenberg 
model [5,17-24,31-42] and for electronic models [25-30]. 

To be specific we consider the Heisenberg antiferromagnet of N spins with quan- 
tum number s = 1/2 in a magnetic field h 

H = J2 J iA ■ *j - hSZ (!) 
and the Hubbard model of N lattice sites 

<*.J> i i,<r=T,i 

<r=t,l 

In (H|) and @ the first sum runs over all neighboring sites on the lattice under 
consideration, Jy > is the antiferromagnetic isotropic Heisenberg exchange inter- 
action between the sites i and j, and S* = J2i s f 1S the z-component of the total 
spin. In ^ tij > is the hopping matrix element between the nearest-neighbor 
sites i and j, U > is the on-site Coulomb repulsion, /j is the chemical potential, 
and n i a = c[ a c i a . For electronic models the chemical potential fi plays the role 
of the magnetic field h. While for the Heisenberg antiferromagnet h controls the 
magnetization M — S z , fi controls the average number of electrons in the system 
for the Hubbard model. 

In what follows we first consider the frustrated quantum Heisenberg antiferro- 
magnet on the sawtooth chain (Fig. [T|) and discuss some generic properties of the 
model which are caused by the localized magnon states. In particular, we consider 
the magnetization process, calculate the ground-state degeneracy of the localized 
eigenstates leading to a finite residual entropy and discuss the low-tcmpcrature ther- 
modynamics for magnetic fields in the vicinity of the saturation field (Seed]). Then 
we illustrate the application of the concepts elaborated for the spin system to the 
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Fig. 2. Left: One-magnon dispersion for the spin-1/2 Heisenberg antiferromagnet on the sawtooth 
chain with Ji = 2, 3\ = 1 and h = (cf. Eq. (0). Right: Ground-state magnetization curves 
m(h) = M(h)/M max for the spin-1/2 Heisenberg antiferromagnet on the sawtooth chain for various 
values of J2 and J\ = 1. 

Hubbard model on the sawtooth chain in Sec. [3] Sec. 2] presents a short summary 
of our discussion. 



2. Localized Magnon States in the Heisenberg Antiferromagnet on 
the Sawtooth Chain 

2.1. Flat bands and localized eigenstates 

In this section we illustrate how the localized magnon states emerge for the frus- 
trated quantum Heisenberg antiferromagnet ([TJ) . The fact that S z commutes with 
the Hamiltonian (fT]) permits us to consider the eigenstates separately in each sub- 
space with different values of S z = N/2, N/2—1, .... In the subspace with S z = N/2 
the only eigenstate is the fully polarized ferromagnetic state, |FM) = | . . 

which plays the role of the vacuum state for the magnon excitations. 

In the one-magnon subspace (S z = N/2 — 1) it is simple to calculate the eigen- 
states given by |1«) = E/U c ' Efii e iKJ s 2 " +/ |FM); H\1 K ) = e±(*s)|l«). The two 
one-magnon branches are given by 

Ji + 2 J 2 1 



e+(k) = h — 
±y ' 2 2 



Jl COS K ± 



\J J\ (-1 + cos nf + 2 Jf (1 + cos k) 



(3) 



For J2 = 2Ji, the lower magnon band becomes completely flat, i.e. £-(/«) = e_ = 
h — 4Ji, see the left panel of Fig. O Let us now focus on the case J2 = 2J\. 
A dispersionless band allows one to construct localized excitations given here by 
|llm) = 4j|FM), where l\- = (1/a/6) ( s 2j-i ~ ^ s 2j + s 2j+i) creates a spin excita- 
tion (magnon) localized in a valley (trapping cell) indicated by bold lines in the 
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upper part of Fig. Q] Note that a typical geometrical feature of a lattice leading to 
the possibility to localize eigenstates is a triangular configuration of antiferromag- 
netic bonds, where the triangle is built by one bond of the trapping cell (here a 
valley) and two bonds attached to the trapping cell [18, 19], see Fig. [U 

Let us consider the n-magnon subspace with S z = N/2 — n. In this subspace 
the construction of the eigenstates of the Heisenberg model is, generally, a difficult 
many-body problem. However, for a lattice which supports localized magnon states, 
a state |nlm) consisting of n independent (i.e. isolated) localized magnons is an ex- 
act eigenstate of the Hamiltonian ([1]). Using the /^-operators introduced above these 
states can be written as |nlm) = l\ l\ ■ ■ ■ l\ |FM), where the ii are sufficiently sepa- 
rated lattice sites. For the sawtooth chain all n localized magnons are independent 
(isolated) if they do not occupy neighboring valleys (hard core rule). This constraint 
immediately leads to a maximum number of localized magnons n max = N/4. The 
energy of the n-particle state |nlm) is 

N 

E n \m = E FM - —h + n(h-4Jx), (4) 

i.e. at h = h% = 4J\ all localized magnon states are degenerate. It is important 
to note that the localized magnon states are the lowest eigenstates in all sectors 
of S z = N/2 - 1, N/2 - 2, ... , N/2 - n max [17, 31]. Hence these states become 
ground states in an appropriate magnetic field. Furthermore it can be shown that 
all localized magnon states are linearly independent for the sawtooth chain [38] and 
that the localized magnon states present the complete manifold of ground states in 
all relevant sectors of S z [21,22,24]. 

In the following sections we will discuss how the localized eigenstates influence 
the physical properties of frustrated lattices. 

2.2. Plateaus and jumps in the magnetization curve 

First we consider the relevance of the localized magnon states for the magnetization 
process. For the calculation of the magnetization M = S z at T = it is sufficient to 
find the lowest energy levels E(M) in the subspaces with different M = N/2, N/2 — 
1, . . . for h = 0. The energy in the presence of an external magnetic field h is given by 
E(M, h) — E(M) — hM, where the magnetization M should acquire a value which 
minimizes E(M, h). Hence M can be determined from the equation dE(M) / dM = h 
which finally gives the magnetization curve m(h) where m = M/M max , M max = 
N/2. For a classical non- frustrated Heisenberg antiferromagnet one typically finds 
a parabolic relation E{M) oc M 2 resulting in a straight-line behavior M oc h. 
Often quantum fluctuations lead only to small deviations from a linear M — h 
relation, see, e.g., Refs. [4,5,43]. However, in the presence of frustration and quantum 
fluctuations more exotic magnetization curves, e.g., curves with plateaus, can be 
observed [4,5,43]. Another spectacular feature observed in magnetization curves of 
frustrated quantum spin systems consists in discontinuous jumps related to a linear 
relation E oc M [5,17-19,24,44]. As discussed in the previous section we find such 
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a linear E — M relation for the sawtooth Heisenberg antiferromagnet with J2 = 
2Ji for values of the magnetization for which the lowest eigenstates are localized 
states, see Eq. (j4j. This leads to a magnetization jump from m = 1/2 directly to 
saturation m = 1 at the saturation field hi = 4Ji, see the right panel of Fig. [2] In 
addition there is wide plateau preceding the jump. This plateau state represents a 
regular pattern of alternately occupied and empty valleys and is two-fold degenerate. 
Magnetization curves with a jump to saturation for other lattices can be found, e.g., 
in Refs. [17-19,33,37,40,43-45]. We emphasize that the jump is macroscopic, and 
that there is no finite-size effect. Furthermore we mention that a jump to saturation 
can be found also for the sawtooth Heisenberg antiferromagnet with higher spin 
quantum number s > 1/2. However, the height of the jump decreases with 1/s, i.e. 
the jump is a true quantum effect and disappears in the classical limit s — > 00. 

Finally, let us discuss deviations from the ideal parameter constellation J2 = 
2 Ji for which the localized magnon states are true eigenstates. The right panel of 
Fig. [5] shows that small deviations (e.g., J2 = 1.9Ji, J2 = 2.1Ji) do not change the 
magnetization curve drastically, whereas the model with uniform bonds J2 = J\ 
exhibits a qualitatively different m(h) behavior. 

2.3. Ground-state residual entropy and low-temperature 
thermodynamics 

It has been shown above that the energy of the n-magnon state in a magnetic field 
is .Efm — Nh/2 + n(h — 4Ji), cf. Eq. (j4|). Obviously, for n < n max this energy 
level is highly degenerate, since there are many ways to place n independent lo- 
calized magnons on a lattice. The degeneracy further increases at the saturation 
field hi = 4Ji, since the energies of the states with different numbers of localized 
magnons n = 0, 1, . . . , n max become equal, namely Eym — Nh±/2. We denote this 
degeneracy at h — h\ by W. Since all localized magnon states are linearly inde- 
pendent [38], they span a highly degenerate ground-state manifold at ft, = h\. The 
degree of degeneracy can be calculated by taking into account the hard-core rule 
(simultaneous occupation of neighboring valleys by localized magnons is forbidden) . 
The remaining counting problem can be solved by mapping the localized magnon 
problem on the sawtooth chain with N sites onto a hard-dimer problem (simulta- 
neous occupation of neighboring sites by dimers is forbidden) on a simple linear 
chain with M = N/2 sites, see the lower part of Fig. [T]and also Refs. [20-22,24]. 
Taking the number of hard-dimer distributions from the literature [46] we can use 
this mapping to find the ground-state degeneracy at the saturation field W. For 
N — * 00 one finds W = ((1 + V / 5)/2) w exp (0.48127V) leading to a finite residual 
entropy of S/k B N = (1/2) In ((1 + VE)/2) w 0.2406 for the sawtooth chain with 
J 2 = 2Ji &th = hi [20-22]. 

In addition, we can use the correspondence between the localized magnon states 
and the spatial configurations of hard dimers to calculate the contribution of the 
localized magnon states to the thermodynamic quantities following the lines given, 
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e.g., in Refs. [46,47]. This contribution may dominate the low-temperature thermo- 
dynamics and therefore we may find predictions for the low-temperature behavior 
of the magnetic quantities in the vicinity of the saturation field hi- The contribution 
of the localized states to the partition function of the spin model can be written as 



Zhn(T, ft, N) = exp 



Eh 



FM 



ftf 



exp 



k B T 
E FM *•* 



z2 9n( 



n) exp 



n=0 



hi 



k B T 



k B T 



S(1>,A0; A/" 



N 
~2' 



(5) 



Here <7Ar(n) is the degeneracy of the ground state of the spin model with ./V sites 
in the sector with n localized magnons, i.e. with M — S z = N/2 — n. In the hard- 
dimer description gN{n) corresponds to the canonical partition function Z(n,J\f) 
of the classical hard-dimer model. h\ — h — pt is the chemical potential of the 
hard dimers and S(T,//,jV) (or E(z,J\f), z = exp (fi/kBT)) is the grand-canonical 
partition function of the one-dimensional hard-dimer lattice gas given by 



1 



exp a:, x 



knT' 



(6) 



Formula (J5J) describes the low-temperature thermodynamics of the spin model 
near the saturation field accurately, i.e. Z(T,h,N) w Z\ m (T,h, N), because of 
the huge degeneracy of the ground state at h = hi (note that there are no 
other ground states apart from the considered localized-magnon states in the cor- 
responding sectors of S z ). We mention that similar considerations are possible 
for other frustrated lattices [21-24,39,41,42]. The contribution of the localized 
magnon states to the Helmholtz free energy F of the spin model is given by 
F lm (T, h, N)/N = E FM /N - ft/2 - k B T\nZ(z,Af)/N. The entropy S, the specific 
heat C, the magnetization M and the susceptibility x follow from F\ m (T, ft, N) 
according to usual relations S hn (T,h,N) = -dF\ m (T, ft, N)/dT, C lm (T,h,N) = 
TdSi m (T,h,N)/dT, M lro (T,ft,iV) = N/2 - (n) = N/2 - k B TdlnE{T, ii,M)/dp, 
X\m(T, ft, N) = dM lm (T, ft, N) /dh. In the limit N -> oo this leads to [20, 23, 24] 



S im (T,h,N) = 1 
k B N 2 

C lm (T,h,N) 



k B N 
M\ m (T, ft, N) 



N 
2 



k B T X i m (T, ft, N) 
N 




1 x exp x 
16 Q+expz) 1 ' 



= 1 - 



+ expx 



exp a; 



16 (i+exp: 



k B T 



hi - h 
k B T 



(7) 
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(h - h^/kjjT (h - hiVkfiT 

Fig. 3. (Color online) The specific heat (left) and the entropy (right) in dependence on x = 
(h± — h)/k^T for the spin-1/2 Heisenberg antiferromagnet on the sawtooth chain with N = 20, 
Ji = 1, J2 = 2 and kgT = 0.01, 0.1, 0.3 in comparison with the one-dimensional hard-dimer (HD) 
gas, Eq. (0. 

The thermodynamic quantities depend on T and h via the universal parameter 
x = (hi — h)/k^T only. Corresponding formulas for finite systems can be found 
using S(T, /j,, AT) from Eq. ([6]) in combination with the relation between F\ m (T, h, N) 
and S(T,fj,,Af) given above. 

Fig. [3] shows a comparison of the entropy and the specific heat of the spin model 
in dependence on the universal parameter x — (hi — h)/kBT with the hard-dimer 
formulas. In addition we show the specific heat as an important measurable quantity 
in dependence on the temperature for magnetic fields slightly above and below the 
saturation field in Fig. 0J 

We emphasize here some prominent features: an extra low-temperature peak in 
the dependence C vs. T for fields slightly below or slightly above h\ (Figs. [3] and [4} 
and an enhanced entropy at h-y at low temperatures (Fig. [3j . Note that C in Eq. {2} 
is zero at x = and consequently there is no extra peak in C(T) for h = hi, see also 
Fig. [3] (left). Furthermore from Figs. [3] and |4] it becomes evident that the hard-dimer 
description works excellently for temperatures up to 10% of the exchange coupling 
and reproduces qualitatively the characteristic features of the spin model for higher 
temperatures up to about 0.3Ji. 

Similar as for the magnetization curve we consider now the influence of de- 
viations from the ideal parameter constellation J2 = 2Ji (for which the localized 
magnon states are true exact eigenstates) on thermodynamic quantities. From Fig. 2] 
it is obvious that only large deviations suppress the extra low-temperature peak in 
C(T). This behavior can be explained by inspection of the low-energy spectrum. 
For small deviations the energy is only slightly changed and the originally highly 
degenerate ground-state manifold becomes quasi-degenerate. As a result the <5-peak 
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h=0.95h, 




Fig. 4. (Color online) Temperature dependence of the specific heat for the spin-1/2 Heisenberg 
antifcrromagnet on the sawtooth chain with N = 20, Ji = 1 and various J2 for magnetic fields 
slightly above (h = X.05hi, left) and below (h = 0.95/ii, right) the saturation field (hi = 4, 4.2, 
3.9 and 3 for J2 = 2, 2.1, 1.9 and 1, respectively). For comparison we show the results for the 
one-dimensional hard-dimer (HD) gas with Af = 10 sites. 

present in the low-energy density of states for J 2 = 2 J\ is broadened but there is 
still a well pronounced maximum in the density of states leading to the extra low-T 
peak in C{T). 

Let us very briefly discuss an aspect of the localized magnon scenario which 
might have some relevance for a possible application of highly frustrated magnets. 
Due to the huge degeneracy of the localized magnon states and the resulting resid- 
ual entropy at ft- = h\ there is a well pronounced low-temperature peak in the 
entropy S versus field h curve, see Fig. [3] (right). It has been pointed out first by 
Zhitomirsky [48] considering the classical kagome Heisenberg antiferromagnet that 
such a degeneracy leads to an enhanced magnetocaloric effect. Later on this point 
has been discussed for quantum spin systems, e.g., in Refs. [20,23,24,49]. 



3. Hubbard Electrons on the Sawtooth Chain 

3.1. Flat one-electron band and localized electron eigenstates 

We consider now the Hubbard model ^ on a sawtooth chain. The specific Hamil- 
tonian reads 



5=0 <r=T4 



ticl^ a c 2j+ 2,a + h \c\j,a C V+l,e + c 2j + 1,<t c 2j+2,<t) + h.C. 
2 1 

+ti(naj,a+ri2j+i,< r ) + U ^ (n 2j ,in 2j ,i + n 2j+ i^n 2j+1A ) , (8) 

3=0 
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where t\ > and t% > are the hopping integrals along the base line and the zig-zag 
path, respectively (see the upper part of Fig. [I]), and U > is the on-site Coulomb 
repulsion. The sawtooth-chain Hubbard model has attracted much attention since 
the 1990s [50]. Here we focus on a special aspect, namely the existence of localized 
ground states and their consequences for the low-temperature physics of the model. 
On the one-particle level the description of the electron system is the same as of 
the XY spin system [26-28]. The one-electron dispersion reads 

e± (k) = (i + t\ cos k ± \Jt\ cos 2 k + 2t\ (1 + cos k) . (9) 

Thus, if ti = V%ti the lowest single electron energy becomes e_ = \i — 2t\, i.e. it is 
completely flat. Similar as for the Heisenberg model we can construct N localized 
one-electron ground states, given by CT |0), l\ ] a = (l/2)(c^-_ 1 a - v^cJy ,*+ c lj+i,*) 
(i.e. the electron is localized in any of the N/2 valleys labeled by the index 2j and 
having either spin up or spin down) with energy e_ = — 2ti + \x. Note that the 
indices of the V and c' operators correspond to the lattice sites as illustrated in 
Fig.ffl 

The Hubbard repulsion becomes relevant in the two-electron subspace. Obvi- 
ously, a two-particle ground state can be constructed by two independent localized 
electrons with arbitrary spin trapped on two valleys which do not touch each other. 
However, in contrast to the Heisenberg model there is no 'hard-core rule', i.e. there 
are further two-particle ground states with two electrons trapped on two neighboring 
valleys, e.g., with indices 2j and 2j + 2. The energy of the corresponding eigenstates 
4j V2j+2 T |0) an d 4j i 4.7+2 1 10) is also independent of U, since both electrons have 
the same spin and therefore the Pauli principle forbids the simultaneous occupation 
of the site 2j + 1 belonging to both valleys. In addition, a straightforward direct 
calculation shows that for two electrons having different spin the linear combination 

4 J -, T 4 J - + 2,jo)+4 J . ; 4 j+2iT |o) ) (io) 

is also a ground state in the two-electron subspace with an energy independent of 
U. This can be seen also by using the SU(2) symmetry of the Hubbard Hamiltonian: 
the state pH|) and the states l\ - |4j+2 t^) an< ^ 4j l4j+2 f° rm a triplet, i.e. ([TO)) 
can be obtained by acting with the total spin lowering operator S~ = c\ , c^j on 

the state ^4j+2 \\^)- ^ course, all states belonging to one triplet have the same 
energy 2(/i - 2t x ). 

We can generalize this procedure to construct the ground states in the subspaces 
with n = 3, . . . , N/2 electrons 

I^WLt-'-'LtI ) ! H\<pi)=n(-2t 1 +n)\<pl). (11) 

They are all degenerate for \i = fio = 2t\ and do not feel U. Evidently, they are 
fully polarized 

S°Wi) = lWl) ; &\ V l) = \ (2 (12) 
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Fig. 5. (Color online) Left: Hole concentration n^/N = 2 — n/N versus chemical potential fj, 
for f2 = \/2ti (localized-electron regime), ti = 0.95\/2ii and t2 = ti for a finite sawtooth chain 
of N = 16 sites (periodic boundary conditions) and U — * oo, t\ = 1. Right: Charge gap Afj, = 
E(N/2 + 1) - 2E(N/2) + E(N/2 - 1) at quarter filling versus U for f 2 = \/2t\, t% = 1 and 
JV= 12,16,20. 

Again the application of S 1- yields new eigenstates with the same energy 
and the same S 2 , but with S z {S~) k \ipl) = (n/2 - fc)(S'-) fc |^). Note that 
4j T " * ^2i fe i ' ' ' 4i! t|0): where i±, . . . ,ik, ■ ■ ■ ,i n denote n contiguous valleys, is not 
an eigenstate. Since there is no hard-core rule the maximum filling with localized 
electrons is n max = N/2, i.e. it is twice as large as for localized magnons. 

In the next step we use the fully polarized n-electron states \cpl) to construct 
the complete set of ground states for < n < N/2. The \<pD can be grouped 
into two classes, namely in one-cluster states and in multi-cluster states. While for 
the one-cluster states the electrons occupy a cluster of contiguous valleys, for a 
multi-cluster state the electrons occupy two or more clusters, where each cluster is 
built by contiguous valleys and the different clusters are separated by one or more 
empty valleys. The key observation is that further ground states can be constructed 
by application of a certain cluster spin flip operator S~ lust — 5Z igclust c\ ,Ci^ on a 
multi-cluster n-electron ground state | The resulting new states are not fully 
polarized and complete the set of ground states in each sector n [28]. 

3.2. Hole concentration in dependence on the chemical potential 

In correspondence to the m{h) curve of the Heisenberg model we consider now the 
hole concentration n-h/N = 2 — n/N in dependence on the chemical potential u 
(Fig. [5] (left)). Like for spin systems, see Sec. 12.21 the main characteristics for the 
system with localized eigenstates (i.e. for ti = v^i) are a size-independent jump 
of rih/N from 3/2 to 2 and a plateau at n^/N = 3/2. This plateau determines the 
range of validity of the localized-electron picture at T = 0. The right panel of Fig. [5] 
presents the plateau width, i.e. the size of the charge gap, versus U for N = 12, 
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16, and 20. One observes that there is almost no finite-size dependence. Since the 
charge gap is zero for U = and increases with U we conclude that its appearance 
is due to the on-site repulsion. Small deviations from the ideal parameter values 
t-2 = \f2t\ do not change the n^/N versus u curve substantially, as illustrated for 
the case t% — 0.95\/2ii in Fig. [5] (left), whereas for the model with uniform hopping 
integrals <2 — ti the charge gap is significantly smaller and there is no indication of 
a jump from the plateau at n^/N = 3/2 to n^/N = 2. 

3.3. Ground-state residual entropy and low-temperature 
thermodynamics 

The localized-electron states are linearly independent, which is connected with the 
fact (as in the case of spin systems, see Ref. [38]) that the middle site is unique to 
each valley. Therefore all these highly degenerate states contribute to the partition 
function. Now the question arises whether the ground state degeneracy can be 
calculated. Due to the different statistics of Hubbard electrons and Heisenberg spins 
there are some differences in the construction rules of localized eigenstates (e.g., the 
occupation of neighboring valleys is forbidden for spins but allowed for electrons, 
see above). Hence it is not surprising that the ground state degeneracy fl'jv(n) for 
n electrons on the TV-site sawtooth chain does not coincide with the one for the 
Heisenberg sawtooth chain (which was equal to the canonical partition functions 
of n hard dimers on a simple chain of AT = N/2 sites, see Sec. 12. 3|) . Nevertheless, 
S'jv(n), n = 0, 1, 2, ... , N/2 for the Hubbard sawtooth chain can also be found by 
a mapping of the localized-electron degrees of freedom onto the one-dimensional 
hard-dimer problem. However, this mapping is more intricate and hard dimers have 
to be considered on a simple chain of ./V sites (instead of N/2 sites as for Heisenberg 
spins), for details see Ref. [28]. One finds <?at(?t.) = Z(n, N) for n = 0, 1, . . . , N/2 — 1 
and g N (N/2) = N/2 + 1 = Z(N/2,N) + N/2 - 1 where Z(n, N) is the canonical 
partition function of the classical one-dimensional hard-dimer model [46,47]. As for 
spin systems we can calculate the contribution of localized electron states to the 
partition function by using this mapping. Again we can present analytical formulas 
for the low-temperature thermodynamic quantities for a non-trivial quantum many- 
body problem. The grand-canonical partition function E of the electron system for 
a chemical potential fi in the vicinity of /iq = 2ti takes the form 



In the thermodynamic limit N — > oo only the largest eigenvalue Ai of the transfer 
matrix survives and, using the definitions S(T,p,N) = ksd{T lnS(T, /i, iV)) /dT, 
C(T,n,N) = TdS(T,fi,N)/dT, we obtain the following results for the thermody- 



S(T>,A0 = Af+A^ + A : 



N 
•3 ' 




(13) 



I 



l 
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Fig. 6. (Color online) Grand-canonical specific heat per site C(T, fi, N)/k^N vs. temperature for 
the sawtooth Hubbard chain of TV = 12 sites for two values of /i, U = oo and t\ = 1, t% = v2, 
0.95^ and 1 (symbols). Note that /i = 2 for t 2 = V2, 0.95^ and 1. For comparison we show 
the hard-dimer data for TV = 12 (solid line) which follows from Eq. (1 1 3 I t and for TV = oo (dashed 
line, Eq. JT4}). Note that for fi = 1.02/^o the hard-dimer data for TV = 12 and TV = oo practically 
coincide. 

namics of one-dimensional hard dimers (see also [28]) 
S(T )fi ,N) _ 




(14) 

4 V i + ex P£ 

which are quite similar to the corresponding expressions for Heisenberg spins, see 
Eq. ([7])- Again we have a finite residual entropy S/k B N = ln((l + \/5)/2) « 0.4812, 
which is twice as large as for the Heisenberg model. 

Results for the low-temperature grand-canonical specific heat are shown in Fig. [6] 
for two values of the chemical potential slightly above and below /i . Similar as for 
the spin system we see (i) that the hard-dimer model, Eqs. (fTo 7 ]) and (fl4|) . yields a 
good description of the electronic model at low temperatures and (ii) that there is 
an extra low-temperature maximum in the grand-canonical specific heat due to the 
manifold of localized electron ground states. Again this additional low-temperature 
maximum in C (T) disappears at \x — \xq as can be read off from Eq. (| 14[) (note that 
C(x = 0) = 0). 

At the end of this section we would like to mention a relation to the so-called 



I 



I 



November 13, 2008 10:29 WSPC/INSTRUCTION FILE paper richter 



14 J. Richter, O. Derzhko, A. Honecker 

flat-band ferromagnetism in the Hubbard model found by Mielke and Tasaki in 
the early 1990s [25]. In particular, the ground states belonging to the plateau at 
n = N/2, see Sec. 13.21 are fully polarized ferromagnetic states. For further details 
of flat-band ferromagnetism in the sawtooth-chain Hubbard model the interested 
reader is referred to the original papers of Tasaki [25] but also to Ref. [28] . 

4. Summary 

To summarize, we have illustrated some basic concepts of localized eigenstates 
in correlated systems on highly frustrated lattices and their effect on the low- 
temperature thermodynamics. As a rule non- interacting electrons or magnons on 
a lattice are delocalized, i.e. are described by a wave function distributed over the 
whole lattice. Electrons or magnons may become localized due to randomness or 
after switching on interactions. As we have discussed on this paper, a frustrating 
lattice topology may lead to another mechanism for localization. Localized states 
may survive in the presence of interactions and under certain conditions they can 
determine the properties of the system at low temperatures. 
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